function [ circularity ] = calcCircularity( vecPhi, vecR )
vecX = vecR.*cos(vecPhi);
vecY = vecR.*sin(vecPhi);

geom = polygeom( vecX, vecY );
%   GEOM = [ area   X_cen  Y_cen  perimeter ]
area = geom(1);
perimeter = geom(4);
circularity = 2*sqrt(pi*area)/perimeter;
end

